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Consider a quantum mechanical system with an unbounded 
hamiltonian. For such a system the entropic index q of non- 
extensive thermodynamics has an upperbound q c > 1 beyond 
which the formalism becomes meaningless. The expression 
l/(q c — 1) is the dimension of the state space (i.e. the man- 
ifold of density matrices) in the context of non-commutative 
geometry. For q = q c an ultraviolet cutoff E < E u ^ is needed 
to guarantee the existence of the equilibrium density matrix. 
Duality between q > 1 and q < 1-statistics via a q — » 1/q 
transformation is established. It leads to an overall picture 
in which the meaning of both q > 1-statistics and q < 1- 
statistics is clarified. The hydrogen atom is considered as an 
example. 
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The aim of this paper is to contribute to the physi- 
cal understanding of nonextensive thermodynamics. It is 
obvious from the many applications (see e.g. |IJ) that the 
formalism is physically relevant. Nevertheless the overall 
picture is still unclear. The open question treated below 
is wether physical systems should be described by an en- 
tropic index q which is less than 1 or larger than 1. The 
q > 1-situation is easier to handle. On the other hand 
there are arguments in favor of q < 1. Both views are 
reconciled by the q — > 1/q-duality which was proposed 
in g]. By working out this duality it becomes feasible 
to attach a physical meaning to both q < 1 and q > 1 
formalisms. 

In H the opinion is expressed that each hamiltonian 
has an intrinsic g-value, and that systems with short 
range interactions all belong to the universality class of 
q = 1-systems. Such an intrinsic value, denoted q c , is 
proposed here. At first sight it does not seem to agree 
with the statement about the universality class of sys- 
tems with short range interactions. However, at the end 
of the paper we will see that a hamiltonian with long 
range potential, like that of the hydrogen atom, does in- 
deed determine a q c -valuc different from 1. 

A side effect of the present work is a link with noncom- 
mutative geometry as formulated in || . This mathemat- 
ical theory generalizes concepts of differential geometry 
to quantum mechanics and could be essential in recon- 
ciling quantum mechanics with theories of gravity. An 
introductory text is B . An essential result of this theory 



is that in a non-commutative context the geometry of 
the manifold of states is fully determined by an operator 
which in in the case of relativistic quantum mechanics is 
the Dirac operator. This correspondence is used here to 
interprete the expression l/(g c — 1) as a geometric dimen- 
sion At this moment it is not clear wether this correspon- 
dence will lead to new insights. It is however not excluded 
that non-extensivity plays an important role in the quan- 
tum mechanics of the microcosmos, as is also suggested 
by previous links between non-extensive thermodynamics 
and quantum groups jq], resp. q-deformations AS97. 
The Tsallis entropy is given by H 



Sg(p) = fcj 



1 - Tr p q 



(1) 



It depends on the entropic parameter q ^ 1 (we assume 
throughout the paper that q > 0) and on the density 
matrix p (p > and Tip — 1). In the limit q — > 1 the 
usual Shannon entropy is recovered. Minimization of the 
free energy 



Ta=U, 
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determines the equilibrium density matrix at tempera- 
ture T. For the energy several proposals have been 
made, the latest of which is 101 



q Kl ' Trp9 



(3) 



The hamiltonian H has a discrete spectrum bounded 
from below. The inverse temperature will be denoted 
(3 = l/ksT as usual. 

Different situations occur wether q > 1 or q < 1. Let 
us first assume that q > 1 . Because of the duality, which 
will be discussed later on, the equilibrium density matrix 
in the q > 1-case is denoted a instead of p. For the same 
reason the hamiltonian is denoted K instead of H . The 
equilibrium density matrix is of the form 



K 



Tro-s 



1/(1-9) 



(4) 



with 



C = Tr 



13(1 -q) 



Trcr9 



1 1/(1-9) 
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Note that Eq. (Q) is of the implicit type. One can prove 
that it has a unique solution, provided that the operator 
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K 1 /( 1 -q) i s trace class (see 0). Let us consider the latter 
condition in more detail. Let (ip n )n be a basis in which K 
is diagonal, with eigenvalues e„: Kip n = e n ip n . Assume 
that the eigenvalues are isolated, with finite multiplicity, 
and ordered increasingly. Let q c be the lower bound of 
q > 1 for which the sequence n 1 ~ q e n is bounded. Then 
obviously the eigenvalues of K diverge as 

e„ ~ n^ 1 (6) 

Note that the dimension of the non-commutative state 
space equals l/(q c — 1). See e.g. [Q, section 5.6. Let us 
consider some examples. If K were bounded (which is 
assumed not to be the case) then q c would be 1. For 
the harmonic oscillator q c equals 2, which corresponds 
with a geometric dimension equal to 1. For a quantum 
particle in a d-dimensional box with reflecting walls q c = 
1 + 2/d. Hence the dimension of the manifold is d/2, 
which is only half of what is expected. The explanation 
is that the laplacean is a second order operator. The 
correct dimension is obtained when the Dirac operator 
or the square root of minus the laplacean is taken as the 
operator determining the geometry of the state space. 

The critical entropic index q c is used as follows. If 
q < q c then IfVCi-?) j s trace class while if q > q c it is 
certainly not. Hence q c is the upperbound of the entropic 
index for which the nonextensive thermodynamics can 
still be formulated (assuming energy has the form (|J)). 

If q = q c then the trace in (^|) diverges generically. 
Finite values can be obtained by introducing a cutoff in 
energy, i.e. by fixing a large number N, and restricting 
the system to the finite Hilbert space spanned by the 
eigenvectors ip n ,n = 0, 1, • • • , N— 1. Note that the cutoff 
in energy E < E u b is related to N by 



E uh ~ 



(7) 



In the limit of large N the trace of (||) can be estimated 
in terms of the Dixmier trace Tr w , which is defined by 



Tr ^.A = lim 

N->ex 



1 



N-l 



In TV 



^<^ n |A|V>„) for all A (8) 



n=0 



The expectation value of any observable A is then given 
by 



Tr al 4- ft' 1 /(i-9c) a 
(A)= lim Tr N a N A= " < * (Q) 

Tr u al + K 1/{1 ~ q ^ 



with a an arbitrary constant such that al 4- K > 0. Re- 
markable is that (A) does not depend on f3. Its value is 
determined by the high energy tail of the spectrum of K. 
In the context of noncommutative geometry (||) is called 
the integral of A. 

The g-averages are given by 



(A) q =Trp N A = irfT ^ 



Tro$ 



(10) 



These are the physical expectation values. They depend 
on temperature and cutoff via a rescaled inverse temper- 
ature 



A straightforward calculation shows that 

^ Tr (/(/?«) + K N )"^ 1 -^A 
\ h- Tr (/(/?*) +i4-)9c/(i-9e) 

with / the inverse of the function g given by 



9(a) 



' Tr (a + K) 



(Tr w (a + i0 1/(1 - ?c) ) 
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(11) 
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Let us now consider the case < q < 1. The equilib- 
rium density matrix is the solution of the implicit equa- 
tion 



1 



1-/3(1 -q) 



H-U q H (p)l 
Trpi 



1/(1-9) 



(14) 



with [A], = A on the subspace on which A is positive, 
and = on the subspace on which A is negative, and 
with 



C = Tr 



P(l -q) 



H-U q H (p)l 
Tvpi 



1/(1-9) 



(15) 



The proof that under certain conditions these equations 
determine a unique equilibrium density matrix is com- 
plicated by the presence of the high energy cutoff — see 
H for a full discussion; part of the problem is discussed 
below. 

In ||lo|| it was argued that a system in thermal contact 
with a finite heath bath is described by Tsallis thermo- 
dynamics with entropic index < q < 1. In this con- 
text the appearance of an energy cutoff is very natural 
since the total energy of system plus heath bath is finite. 
This bound on the total energy belongs probably to the 
essence of nonextensive systems. 

A q ~ 1/g-duality has been studied in M. It is based 
on the observation that the relation p = a q / Tra q can be 
inverted to a — p q / Tr p q with q' = 1/q. This duality is 
worked out below by introduction of the effective Hamil- 
tonian _ff c ff of p| . Write the energy functional ([|) in a 
1-homogeneous form 



Wf(p) = (Trp) 



Tr p q H 
Tip q 



(16) 



Then the effective Hamiltonian iJ c ff is defined by varia- 
tion of the energy functional 



H cfi (p) = j-U?{p) 



(17) 
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(the variation is restricted to the subspace on which p is 
strictly positive). It is the generator of the time evolution 
(see Jllj ) . For an equilibrium state [p, H ] = holds, in 
which case it is easy to obtain (putting Trp = 1) 

H oS (p) =U^(p)l + ^- qP ^ 1 (H-U q H (p)l) (18) 

Note also that 

Tr oH = TvpH cS {p) = U» \p) (19) 

because of the 1-homogeneity of (p!q). A short calculation 
now gives 



1-/3(1-9) 
1 



Trp? 



AV g Tr ct / 9 



(20) 



with A = TrcrVv/Trp". The lhs of @ is proportional 
to p 1-9 on the subspace on which it is strictly positive. 
Hence from a ~ p 9 follows 



/3(1 



1 H eS (p)-U?(p)l 



1/(1-1/9) 



with £ the normalization factor. It turns out that 
£ = C 9 Trp 9 



(21) 



(22) 



There is no cutoff in ( |2l| ) because H Q g(p) is zero on the 
subspace where the cutoff is active. ( pl| ) shows that <r 
is the equilibrium density matrix corresponding to the 
entropic index 1/q for a system with hamiltonian K — 
\H c ff (p) . One concludes that a nice duality between q < 
1 and q > 1-statistics is obtained if not only q is mapped 
onto 1/q and p is mapped onto a, but also H is mapped 
onto K — XH c g(p). 

From this duality the following statements can be de- 
duced. 

• q > 1-statistics is another way of looking at q < 1- 
statistics. The population of energy levels as de- 
scribed by the q < 1-density matrix p drops faster 
than exponentially with increasing energy, as it is 
expected in a non-extensive system in thermal con- 
tact with a finite heath bath. The density matrix <j 
of q > 1-statistics is a renormalized density matrix 
in which the population of high energy levels has 
been increased in an artificial way. 

• Not every system described with q > 1-statistics 
can be obtained in this way. Indeed, due to the cut- 
off in (|l4|) the effective hamiltonian H & s is in most 
cases defined on a finite dimensional Hilbert space. 
Hence a cutoff should be applied to the q > 1- 
hamiltonian K before considering it as an effective 



hamiltonian. Further complications are that the 
effective hamiltonian depends on temperature, and 
that it is not the result of a simple ultraviolet cutoff 
applied to a fixed hamiltonian. 

• Two types of hamiltonians are involved. The hamil- 
tonian H in the q < 1-description is the operator 
which determines the energy of the system. The op- 
erator K of the q > 1-description corresponds with 
the effective hamiltonian H e ff of the q < 1-case and 
hence is also the generator of the dynamics. This 
observation is in agreement with the use of K in the 
context of noncommutative geometry to formulate 
Tomita-Takesaki-theory (see ||, section 5.7). The 
connection with linear respons theory and the KMS 
condition will be discussed elsewhere (l2) . 

Let us now consider which kind of system at q < 1 cor- 
responds by duality with a system with entropic index q 
satisfying 1 < q < q c . For most hamiltonians the equi- 
librium density matrix of the < q < 1-case will contain 
only a finite number of eigenvalues different from zero. 
This is a consequence of the high-energy cutoff in (0). 
However, there is one exception, when the hamiltonian 
H is bounded and has an accumulation of eigenvalues 
at the upper limit of its spectrum. This is e.g. the 
case with the discrete part of the spectrum of the hydro- 
gen atom. This is precisely the kind of hamiltonian for 
which q < 1-statistics is physically meaningful because 
the usual Gibbs statistics cannot be applied - see p3[ . 
Let f3- be the solution of the equation 



Trp« = /3(l-g)( eoo -Wf(p)) 



(23) 



If (3 > j3- then the energy cutoff is active and only a finite 
number of eigenvalues of the equilibrium density matrix 
differ from zero. At high temperatures, i.e. f3 < /3_, the 
operator in (Q) has a diverging trace and g-statistics 
is meaningless, as it is the case for q = 1. The inverse 
temperature /?_ depends on q. It tends to oo when q — > 1 
and when q — > l/q c , with the critical entropic index q c 
equal to the upperbound of q's for which n q ~ 1 (e QO — e n ) 
is bounded, i.e. 



1 



n 9c 



(24) 



One concludes that the thermodynamic equilibrium of 
q < 1-statistics exists for q in the interval (l/q c , 1) and for 
low enough temperatures so that f3 > (3-. For complete- 
ness, note that for a bounded hamiltonian in an infinite 
dimensional Hilbert space the free energy has no mini- 
mum. Indeed, the energy U q (p) is bounded above, while 
the entropy S q (p) can be made arbitrary large. Hence 
the free energy !F q {p) can be made arbitrary small. The 
density matrix given by ( pT[ ) is only a relative minimum 
of the free energy. However it is the unique maximum of 



3 



the entropy S q (p) for a given value of the energy U q (p) 
-see|. 

Consider the example of the discrete part of the spec- 
trum of the hydrogen atom. It has g c = 5/3 (indeed the 
spectrum consists of eigenvalues —a/n 2 , n = 1, 2, • • •; the 
n-th level is ro 2 -times degenerate; due to the degeneracy, 
the eigenvalues e„ tend to — as —a/n 2 / 3 ; neglect 
the spin of the electron). The corresponding dimension 
equals 3/2 (again, one would expect the dimension to be 
equal to 3; the extra factor 1/2 disappears if the square 
root of minus the hamiltonian is the operator determin- 
ing the geometry). See the figure for a plot of a/3_ as a 
function of q in the domain of existence (3/5, 1). 
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FIG. 1. aft- as a function of q. 

In Jl3| an anomaly is reported at q c — 9/7 instead 
of q c — 5/3. This value was taken from JL4| , which is 
a paper on stellar polytropes. Study of the latter paper 
makes clear that it concerns the same anomaly as the one 
studied here, namely the upper limit for which ^-statistics 
is meaningful. However, it is not evident that q c should 
have the same value in a classical model of gravitation as 
in a quantum mechanical model of Coulomb attraction. 
Further differences between |l3j and the present paper 
can be explained by the use of a different expression for 
the energy functional (^|). 

One concludes that a statistical description of the hy- 
drogen atom is possible in the interval q~ x < q < 1. The 
temperature varies between zero for the ground state and 
the maximal value T_ = l/ksfl- which depends itself on 
q. If the heath bath is very large then q is close to one and 
the abundance of energy destabilizes the system (i.e. T_ 
tends to zero). If q tends to q c then the Tsallis entropy 
diverges and again T_ tends to zero. 
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